Epigenetic variation impacts individual differences in the transcriptional response to influenza infection.

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Characterization of caffeine response regulatory variants in vascular endothelial cells.

Genetic variants in gene regulatory sequences can modify gene expression and mediate the molecular response to environmental stimuli. In addition, genotype-environment interactions (GxE) contribute to complex traits such as cardiovascular disease. Caffeine is the most widely consumed stimulant and is known to produce a vascular response. To investigate GxE for caffeine, we treated vascular endothelial cells with caffeine and used a massively parallel reporter assay to measure allelic effects on gene regulation for over 43,000 genetic variants. We identified 665 variants with allelic effects on gene regulation and 6 variants that regulate the gene expression response to caffeine (GxE, false discovery rate [FDR] < 5%). When overlapping our GxE results with expression quantitative trait loci colocalized with coronary artery disease and hypertension, we dissected their regulatory mechanisms and showed a modulatory role for caffeine. Our results demonstrate that massively parallel reporter assay is a powerful approach to identify and molecularly characterize GxE in the specific context of caffeine consumption.

Deciphering maternal-fetal cross-talk in the human placenta during parturition using single-cell RNA sequencing.

Labor is a complex physiological process requiring a well-orchestrated dialogue between the mother and fetus. However, the cellular contributions and communications that facilitate maternal-fetal cross-talk in labor have not been fully elucidated. Here, single-cell RNA sequencing (scRNA-seq) was applied to decipher maternal-fetal signaling in the human placenta during term labor. First, a single-cell atlas of the human placenta was established, demonstrating that maternal and fetal cell types underwent changes in transcriptomic activity during labor. Cell types most affected by labor were fetal stromal and maternal decidual cells in the chorioamniotic membranes (CAMs) and maternal and fetal myeloid cells in the placenta. Cell-cell interaction analyses showed that CAM and placental cell types participated in labor-driven maternal and fetal signaling, including the collagen, C-X-C motif ligand (CXCL), tumor necrosis factor (TNF), galectin, and interleukin-6 (IL-6) pathways. Integration of scRNA-seq data with publicly available bulk transcriptomic data showed that placenta-derived scRNA-seq signatures could be monitored in the maternal circulation throughout gestation and in labor. Moreover, comparative analysis revealed that placenta-derived signatures in term labor were mirrored by those in spontaneous preterm labor and birth. Furthermore, we demonstrated that early in gestation, labor-specific, placenta-derived signatures could be detected in the circulation of women destined to undergo spontaneous preterm birth, with either intact or prelabor ruptured membranes. Collectively, our findings provide insight into the maternal-fetal cross-talk of human parturition and suggest that placenta-derived single-cell signatures can aid in the development of noninvasive biomarkers for the prediction of preterm birth.

Blockade of IL-6R prevents preterm birth and adverse neonatal outcomes.

BACKGROUND: Preterm birth preceded by spontaneous preterm labour often occurs in the clinical setting of sterile intra-amniotic inflammation (SIAI), a condition that currently lacks treatment. METHODS: Proteomic and scRNA-seq human data were analysed to evaluate the role of IL-6 and IL-1α in SIAI. A C57BL/6 murine model of SIAI-induced preterm birth was developed by the ultrasound-guided intra-amniotic injection of IL-1α. The blockade of IL-6R by using an aIL-6R was tested as prenatal treatment for preterm birth and adverse neonatal outcomes. QUEST-MRI evaluated brain oxidative stress in utero. Targeted transcriptomic profiling assessed maternal, foetal, and neonatal inflammation. Neonatal biometrics and neurodevelopment were tested. The neonatal gut immune-microbiome was evaluated using metagenomic sequencing and immunophenotyping. FINDINGS: IL-6 plays a critical role in the human intra-amniotic inflammatory response, which is associated with elevated concentrations of the alarmin IL-1α. Intra-amniotic injection of IL-1α resembles SIAI, inducing preterm birth (7% vs. 50%, p = 0.03, Fisher's exact test) and neonatal mortality (18% vs. 56%, p = 0.02, Mann-Whitney U-test). QUEST-MRI revealed no foetal brain oxidative stress upon in utero IL-1α exposure (p > 0.05, mixed linear model). Prenatal treatment with aIL-6R abrogated IL-1α-induced preterm birth (50% vs. 7%, p = 0.03, Fisher's exact test) by dampening inflammatory processes associated with the common pathway of labour. Importantly, aIL-6R reduces neonatal mortality (56% vs. 22%, p = 0.03, Mann-Whitney U-test) by crossing from the mother to the amniotic cavity, dampening foetal organ inflammation and improving growth. Beneficial effects of prenatal IL-6R blockade carried over to neonatal life, improving survival, growth, neurodevelopment, and gut immune homeostasis. INTERPRETATION: IL-6R blockade can serve as a strategy to treat SIAI, preventing preterm birth and adverse neonatal outcomes. FUNDING: NICHD/NIH/DHHS, Contract HHSN275201300006C. WSU Perinatal Initiative in Maternal, Perinatal and Child Health.

Single nucleus transcriptomics of ventral midbrain identifies glial activation associated with chronic opioid use disorder.

Dynamic interactions of neurons and glia in the ventral midbrain mediate reward and addiction behavior. We studied gene expression in 212,713 ventral midbrain single nuclei from 95 individuals with history of opioid misuse, and individuals without drug exposure. Chronic exposure to opioids was not associated with change in proportions of glial and neuronal subtypes, however glial transcriptomes were broadly altered, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes. Genes associated with activation of the immune response including interferon, NFkB signaling, and cell motility pathways were upregulated, contrasting with down-regulated expression of synaptic signaling and plasticity genes in ventral midbrain non-dopaminergic neurons. Ventral midbrain transcriptomic reprogramming in the context of chronic opioid exposure included 325 genes that previous genome-wide studies had linked to risk of substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

Dissecting human population variation in single-cell responses to SARS-CoV-2.

Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution.

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Single Nucleus Transcriptomics Reveals Pervasive Glial Activation in Opioid Overdose Cases.

Dynamic interactions of neurons and glia in the ventral midbrain (VM) mediate reward and addiction behavior. We studied gene expression in 212,713 VM single nuclei from 95 human opioid overdose cases and drug-free controls. Chronic exposure to opioids left numerical proportions of VM glial and neuronal subtypes unaltered, while broadly affecting glial transcriptomes, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes, with prominent activation of the immune response including interferon, NFkB signaling, and cell motility pathways, sharply contrasting with down-regulated expression of synaptic signaling and plasticity genes in VM non-dopaminergic neurons. VM transcriptomic reprogramming in the context of opioid exposure and overdose included 325 genes with genetic variation linked to substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

The amniotic fluid proteome changes with term labor and informs biomarker discovery in maternal plasma.

The intra-uterine components of labor, namely, myometrial contractility, cervical ripening, and decidua/membrane activation, have been extensively characterized and involve a local pro-inflammatory milieu of cellular and soluble immune mediators. Targeted profiling has demonstrated that such processes extend to the intra-amniotic space, yet unbiased analyses of the proteome of human amniotic fluid during labor are lacking. Herein, we utilized an aptamer-based platform to characterize 1,310 amniotic fluid proteins and found that the proteome undergoes substantial changes with term labor (251 proteins with differential abundance, q < 0.1, and fold change > 1.25). Proteins with increased abundance in labor are enriched for immune and inflammatory processes, consistent with prior reports of labor-associated changes in the intra-uterine space. By integrating the amniotic fluid proteome with previously generated placental-derived single-cell RNA-seq data, we demonstrated the labor-driven upregulation of signatures corresponding to stromal-3 and decidual cells. We also determined that changes in amniotic fluid protein abundance are reflected in the maternal plasma proteome. Collectively, these findings provide novel insights into the amniotic fluid proteome in term labor and support its potential use as a source of biomarkers to distinguish between true and false labor by using maternal blood samples.

Probabilistic integration of transcriptome-wide association studies and colocalization analysis identifies key molecular pathways of complex traits.

Integrative genetic association methods have shown great promise in post-GWAS (genome-wide association study) analyses, in which one of the most challenging tasks is identifying putative causal genes and uncovering molecular mechanisms of complex traits. Recent studies suggest that prevailing computational approaches, including transcriptome-wide association studies (TWASs) and colocalization analysis, are individually imperfect, but their joint usage can yield robust and powerful inference results. This paper presents INTACT, a computational framework to integrate probabilistic evidence from these distinct types of analyses and implicate putative causal genes. This procedure is flexible and can work with a wide range of existing integrative analysis approaches. It has the unique ability to quantify the uncertainty of implicated genes, enabling rigorous control of false-positive discoveries. Taking advantage of this highly desirable feature, we further propose an efficient algorithm, INTACT-GSE, for gene set enrichment analysis based on the integrated probabilistic evidence. We examine the proposed computational methods and illustrate their improved performance over the existing approaches through simulation studies. We apply the proposed methods to analyze the multi-tissue eQTL data from the GTEx project and eight large-scale complex- and molecular-trait GWAS datasets from multiple consortia and the UK Biobank. Overall, we find that the proposed methods markedly improve the existing putative gene implication methods and are particularly advantageous in evaluating and identifying key gene sets and biological pathways underlying complex traits.

Analysis of transcriptional changes in the immune system associated with pubertal development in a longitudinal cohort of children with asthma.

Puberty is an important developmental period marked by hormonal, metabolic and immune changes. Puberty also marks a shift in sex differences in susceptibility to asthma. Yet, little is known about the gene expression changes in immune cells that occur during pubertal development. Here we assess pubertal development and leukocyte gene expression in a longitudinal cohort of 251 children with asthma. We identify substantial gene expression changes associated with age and pubertal development. Gene expression changes between pre- and post-menarcheal females suggest a shift from predominantly innate to adaptive immunity. We show that genetic effects on gene expression change dynamically during pubertal development. Gene expression changes during puberty are correlated with gene expression changes associated with asthma and may explain sex differences in prevalence. Our results show that molecular data used to study the genetics of early onset diseases should consider pubertal development as an important factor that modifies the transcriptome.

A single-cell atlas of murine reproductive tissues during preterm labor.

Preterm birth, the leading cause of perinatal morbidity and mortality worldwide, frequently results from the syndrome of preterm labor. The best-established causal link to preterm labor is intra-amniotic infection, which involves premature activation of the parturition cascade in the reproductive tissues. Herein, we utilize single-cell RNA sequencing (scRNA-seq) to generate a single-cell atlas of the murine uterus, decidua, and cervix in a model of infection-induced preterm labor. We show that preterm labor affects the transcriptomic profiles of specific immune and non-immune cell subsets. Shared and tissue-specific gene expression signatures are identified among affected cells. Determination of intercellular communications implicates specific cell types in preterm labor-associated signaling pathways across tissues. In silico comparison of murine and human uterine cell-cell interactions reveals conserved signaling pathways implicated in labor. Thus, our scRNA-seq data provide insights into the preterm labor-driven cellular landscape and communications in reproductive tissues.

Preparation of single-cell suspensions from the human placenta.

Single-cell RNA-sequencing (scRNA-seq) allows the characterization of cellular composition and interactions in complex tissues. An essential prerequisite for scRNA-seq is the preparation of high-quality single-cell suspensions. So far, no protocols have been described for preparing such suspensions from the placenta, an essential organ for fetal development and a site of maternal-fetal immune interaction. Here we describe a protocol for the preparation of high-quality single-cell suspensions from human placental tissues-namely, the basal plate, placental villi and chorioamniotic membranes. The protocol outlines the collection of tissues from the placenta, tailored dissociation procedures for each tissue, and the cryopreservation of single-cell suspensions for multiplex sequencing library preparation. The protocol can be performed by a qualified investigator with basic working knowledge of placental structure. Moreover, the single-cell suspensions generated by using this protocol are compatible with droplet-based scRNA-seq technology, such as the 10x Genomics Chromium system. This protocol reliably produces single-cell suspensions from the placental tissues with high yield and viability for scRNA-seq. This protocol takes ~6 h to complete from tissue collection to cryopreservation of single-cell suspensions, and an additional 2 h for thawing of cryopreserved single cells.

Single-Cell Immunobiology of the Maternal-Fetal Interface.

Pregnancy success requires constant dialogue between the mother and developing conceptus. Such crosstalk is facilitated through complex interactions between maternal and fetal cells at distinct tissue sites, collectively termed the "maternal-fetal interface." The emergence of single-cell technologies has enabled a deeper understanding of the unique processes taking place at the maternal-fetal interface as well as the discovery of novel pathways and immune and nonimmune cell types. Single-cell approaches have also been applied to decipher the cellular dynamics throughout pregnancy, in parturition, and in obstetrical syndromes such as recurrent spontaneous abortion, preeclampsia, and preterm labor. Furthermore, single-cell technologies have been used during the recent COVID-19 pandemic to evaluate placental viral cell entry and the impact of SARS-CoV-2 infection on maternal and fetal immunity. In this brief review, we summarize the current knowledge of cellular immunobiology in pregnancy and its complications that has been generated through single-cell investigations of the maternal-fetal interface.

Human Plasma Proteome During Normal Pregnancy.

The human plasma proteome is underexplored despite its potential value for monitoring health and disease. Herein, using a recently developed aptamer-based platform, we profiled 7288 proteins in 528 plasma samples from 91 normal pregnancies (Gene Expression Omnibus identifier GSE206454). The coefficient of variation was <20% for 93% of analytes (median 7%), and a cross-platform correlation for selected key angiogenic and anti-angiogenic proteins was significant. Gestational age was associated with changes in 953 proteins, including highly modulated placenta- and decidua-specific proteins, and they were enriched in biological processes including regulation of growth, angiogenesis, immunity, and inflammation. The abundance of proteins corresponding to RNAs specific to populations of cells previously described by single-cell RNA-Seq analysis of the placenta was highly modulated throughout gestation. Furthermore, machine learning-based prediction of gestational age and of time from sampling to term delivery compared favorably with transcriptomic models (mean absolute error of 2 weeks). These results suggested that the plasma proteome may provide a non-invasive readout of placental cellular dynamics and serve as a blueprint for investigating obstetrical disease.

Deconstructing cold-induced brown adipocyte neogenesis in mice.

Cold exposure triggers neogenesis in classic interscapular brown adipose tissue (iBAT) that involves activation of ß1-adrenergic receptors, proliferation of PDGFRA+ adipose tissue stromal cells (ASCs), and recruitment of immune cells whose phenotypes are presently unknown. Single-cell RNA-sequencing (scRNA-seq) in mice identified three ASC subpopulations that occupied distinct tissue locations. Of these, interstitial ASC1 were found to be direct precursors of new brown adipocytes (BAs). Surprisingly, knockout of ß1-adrenergic receptors in ASCs did not prevent cold-induced neogenesis, whereas pharmacological activation of the ß3-adrenergic receptor on BAs was sufficient, suggesting that signals derived from mature BAs indirectly trigger ASC proliferation and differentiation. In this regard, cold exposure induced the delayed appearance of multiple macrophage and dendritic cell populations whose recruitment strongly correlated with the onset and magnitude of neogenesis across diverse experimental conditions. High-resolution immunofluorescence and single-molecule fluorescence in situ hybridization demonstrated that cold-induced neogenesis involves dynamic interactions between ASC1 and recruited immune cells that occur on the micrometer scale in distinct tissue regions. Our results indicate that neogenesis is not a reflexive response of progenitors to ß-adrenergic signaling, but rather is a complex adaptive response to elevated metabolic demand within brown adipocytes.

Proteoglycans: Systems-Level Insight into Their Expression in Healthy and Diseased Placentas.

Proteoglycan macromolecules play key roles in several physiological processes (e.g., adhesion, proliferation, migration, invasion, angiogenesis, and apoptosis), all of which are important for placentation and healthy pregnancy. However, their precise roles in human reproduction have not been clarified. To fill this gap, herein, we provide an overview of the proteoglycans' expression and role in the placenta, in trophoblast development, and in pregnancy complications (pre-eclampsia, fetal growth restriction), highlighting one of the most important members of this family, syndecan-1 (SDC1). Microarray data analysis showed that of 34 placentally expressed proteoglycans, SDC1 production is markedly the highest in the placenta and that SDC1 is the most upregulated gene during trophoblast differentiation into the syncytiotrophoblast. Furthermore, placental transcriptomic data identified dysregulated proteoglycan genes in pre-eclampsia and in fetal growth restriction, including SDC1, which is supported by the lower concentration of syndecan-1 in maternal blood in these syndromes. Overall, our clinical and in vitro studies, data analyses, and literature search pointed out that proteoglycans, as important components of the placenta, may regulate various stages of placental development and participate in the maintenance of a healthy pregnancy. Moreover, syndecan-1 may serve as a useful marker of syncytialization and a prognostic marker of adverse pregnancy outcomes. Further studies are warranted to explore the role of proteoglycans in healthy and complicated pregnancies, which may help in diagnostic or therapeutic developments.

Analyzing and reconciling colocalization and transcriptome-wide association studies from the perspective of inferential reproducibility.

Transcriptome-wide association studies and colocalization analysis are popular computational approaches for integrating genetic-association data from molecular and complex traits. They show the unique ability to go beyond variant-level genetic-association evidence and implicate critical functional units, e.g., genes, in disease etiology. However, in practice, when the two approaches are applied to the same molecular and complex-trait data, the inference results can be markedly different. This paper systematically investigates the inferential reproducibility between the two approaches through theoretical derivation, numerical experiments, and analyses of four complex trait GWAS and GTEx eQTL data. We identify two classes of inconsistent inference results. We find that the first class of inconsistent results (i.e., genes with strong colocalization but weak transcriptome-wide association study [TWAS] signals) might suggest an interesting biological phenomenon, i.e., horizontal pleiotropy; thus, the two approaches are truly complementary. The inconsistency in the second class (i.e., genes with weak colocalization but strong TWAS signals) can be understood and effectively reconciled. To this end, we propose a computational approach for locus-level colocalization analysis. We demonstrate that the joint TWAS and locus-level colocalization analysis improves specificity and sensitivity for implicating biologically relevant genes.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single cell resolution.

Synthetic glucocorticoids, such as dexamethasone, have been used as treatment for many immune conditions, such as asthma and more recently severe COVID-19. Single cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated Peripheral Blood Mononuclear Cells from 96 African American children. We employed novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we demonstrated that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and demonstrated widespread genetic regulation of the transcriptional dynamics of the gene expression response.

IL-22 Plays a Dual Role in the Amniotic Cavity: Tissue Injury and Host Defense against Microbes in Preterm Labor.

IL-22 is a multifaceted cytokine with both pro- and anti-inflammatory functions that is implicated in multiple pathologies. However, the role of IL-22 in maternal-fetal immunity in late gestation is poorly understood. In this study, we first showed that IL-22+ T cells coexpressing retinoic acid-related orphan receptor γt (ROR-γt) are enriched at the human maternal-fetal interface of women with preterm labor and birth, which was confirmed by in silico analysis of single-cell RNA sequencing data. T cell activation leading to preterm birth in mice was preceded by a surge in IL-22 in the maternal circulation and amniotic cavity; however, systemic administration of IL-22 in mice did not induce adverse perinatal outcomes. Next, using an ex vivo human system, we showed that IL-22 can cross from the choriodecidua to the intra-amniotic space, where its receptors (Il22ra1, Il10rb, and Il22ra2) are highly expressed by murine gestational and fetal tissues in late pregnancy. Importantly, amniotic fluid concentrations of IL-22 were elevated in women with sterile or microbial intra-amniotic inflammation, suggesting a dual role for this cytokine. The intra-amniotic administration of IL-22 alone shortened gestation and caused neonatal death in mice, with the latter outcome involving lung maturation and inflammation. IL-22 plays a role in host response by participating in the intra-amniotic inflammatory milieu preceding Ureaplasma parvum-induced preterm birth in mice, which was rescued by the deficiency of IL-22. Collectively, these data show that IL-22 alone is capable of causing fetal injury leading to neonatal death and can participate in host defense against microbial invasion of the amniotic cavity leading to preterm labor and birth.